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Abstract 

Observations show that there is significant clustering of QSO metal absorption lines 
within the range of velocity dispersion between 200km/sec and 600 km/sec. With a reason- 
able supernova rate, it is shown that high velocity gases driven by SNe and/or strong stellar 
winds could explain the clustered absorptions, provided that QSO metal-line absorbers are 
galactic halos or dwarf galaxies. Rich clusters of galaxies, on the other hand, cannot yield 
the observed clustering of QSO metal absorption lines. 
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I. INTRODUCTION 



The metal absorption lines in QSO spectra (type C in the classification of Weyman, 
Carswell & Smith (1981)) have raised great interests because they might reveal properties 
and evolution of objects associated with galaxies at redshifts z > 1. The properties of these 
metal absorption fines were reviewed by Weyman, Carswell & Smith (1981) and recently by 
York et al. (1991). These absorption lines are commonly associated with halos of galaxies 
(Bahcall & Spitzer 1969), since halos contain metals, have star formation, and are clustered 
fike galaxies. It is also reasonable if considering that the cross section of the absorbers (if 
they associate with galaxies) needed to explain these lines is ~ 10 times the visible extent 
of modern galaxies (Sargent et al. 1979). Three options for "halos" have been proposed: 
1. each absorber corresponds to a single halo (Pettini et al. 1983); 2. a small number of 
gas-rich dwarf galaxies (like the Magellanic clouds) form each halo (York et al. 1986); 3. 
many mini- absorbers (with masses much less than those of dwarf galaxies, such as lO^M©) 
in a halo provide the metal- line absorption (Sargent et al. 1979). 

A significant portion of these metal-line absorption systems are systems that have mul- 
tiple absorption lines with velocity splits > 200 km/sec, and lie in the intermediate position 
between quasars and observers (Pettini et al. 1983; York et al. 1984; Sargent, Boksen- 
berg, & Steidel 1988, hereafter SBS). A two-point correlation calculation of CIV absorption 
lines in QSO spectra shows that there is significant correlation in the range of 200 km/sec 
< Ay < 600 km/sec, where is the velocity split between a pair of absorption lines 
(SBS). In the survey of 55 QSOs by SBS, 111 CIV absorbers are identified with an equiva- 
lent width larger than 0.15A and a relative velocity with respect to the QSO larger than 5000 
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km/sec; among them 18 pairs of absorbers are found in excess of what would be expected 
from a random Poisson distribution of absorbers in the range of 200 km/sec < AV < 600 
km/sec, while no significant correlation is found in the range of 600 km/sec< AV < 1000 
km/sec. In other words, metal absorption lines are found to cluster in the range of velocity 
dispersion between 200 km/sec and 600 km/sec. 

The clustering of absorption lines in the range of 200 km/sec< AV < 600 km/sec 
cannot, at least mainly, be explained by the differential rotation motion within a galaxy, 
which typically yields a AV of 150 km/sec. Attempts have been made to explain these 
clustered absorptions by galaxy-galaxy correlation at z ~ 2, the average redshift of these 
absorbers (SBS). However, it was shown that the amplitude of today's galaxy-galaxy cor- 
relation function extrapolated to 2; ~ 2, or the galaxy-galaxy correlation calculated from 
cold dark matter model, fell short to explain the amplitude of the clustering among these 
absorbers by a factor >3 (SBSQ Salmon & Hogan 1986). Although a high biasing factor 
may improve the fit of the galaxy-galaxy correlation function to the absorber correlation at 

*the expected correlation amplitude calculated by SBS from their eq. (24) seems too high. Using 
all the notations in SBS, it is clear from their eq. (24) that nc{R) < 2itB? ^{rc/rY'^dr. 
Since $ = $o(l + -2)^ and Tc = 5h~^{l + 2;)~^/^Mpc (SBS), the number of correlated galaxies 
per galaxy divided by the mean number of absorbers per unit redshift, ric/N < 2.5 • 5/i^^Mpc ■ 
{5h'^Mpc/6R^Y■^/{dn^/TT^). With 70/i-^kpc, dn = 3000/i-^Mpc and z ~ 2, njN < 

0.003, a factor of 20 smaller than the 0.062 quoted in SBS. If a gravitationally unbound galaxy- 
galaxy correlation is assumed, i.e., tq = 5(1 + z)~^Mpc, the resultant ric/N < 0.01, which still a 
factor of 6 lower than quoted in SBS. 
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200km/sec< < 600km/sec, it may be hard to understand why a significant correlation 
of absorbers doesn't extend to AV > 600 km/sec, given the high biasing and a not-so-steep 
power law galaxy-galaxy correlation function (Peebles 1993) 



where r is the distance between a pair of galaxies and tq is a constant. 

In this letter, we set aside the galaxy-galaxy correlation explanation to these clustered 
absorption systems, and investigate two other alternatives: 1) clusters of galaxies which 
have velocity dispersions as high as 1000 km/sec intercepting quasar lines of sight; 2) high 
velocity gases inside absorbers driven by supernova shocks and/or stellar winds. We will 
show that while the first alternative fails to explain the observed clustering of QSO metal 
absorption lines at 200km/sec < AV < 600km/sec, the second alternative may provide a 
viable explanation. 



If a quasar line of sight strikes two or more galaxies in a rich cluster, they can provide 
clustered absorptions with a velocity spread as high as ~ 1000 km/sec, which is the typical 
velocity dispersion of a cluster (Pettini et al. 1983). However, since rich clusters only have 
a dense core within a 0.2h~^ Mpc radius (Peebles 1993), the chance of a quasar line of sight 
striking multiple members of a single cluster is small. 



Observations show that average rich clusters have a galaxy density distribution of (Lilje 
& Efstathiou 1988) 




(1) 



II. CLUSTERS OF GALAXIES 



n{r) = nb[l + (ro/r)"^'^], for 0.2/i-^Mpc < r < 



20/i-^Mpc, 



(2) 
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where n(r) is the number density of galaxies at a distance r to the cluster center, 
riTa is the background number density of galaxies, Tq 9/i^^(l + z)^^'^Mpc, and h = 
iJo/lOOkm/sec/Mpc with Hq being the Hubble constant. From APM survey nb ~ 
0.014/i3(l + zf Mpc-3 (Loveday et al. 1992). n{r) falls off much faster beyond 20/i~^Mpc 
(Lilje & Efstathiou 1988), and is much smoother inside the cluster core. In our following 
calculations, it is convenient to assume a uniformly distributed core with 

n(r) = ncore < (2i?h)"', for r < 0.2/i-^Mpc, (3) 

where i?h is the typical radius of a galactic halo. The detailed functional form of n(r) in the 
core, however, does not affect the generality of our calculations in any significant manner. 
If each galactic halo provides each individual absorber, to explain the number of metal-line 
absorption systems observed, it was shown that (SBS) 

i?h ~ 70h-\pc. (4) 

Therefore, when a quasar line of sight passes through the cluster with the closest distance 
r_L from the cluster center (as shown in Figure 1), the number of galaxies that the quasar 
line of sight would intercept is 

271" /o " ^ R^n{r)dr', for r± > Tc] 

^('■^) = { rrrr (5) 
27r Rln,,r.dr' + l-n / Rln{r)dr' , for < r,. 



In eq. (|^), r = yr'^ + r^^. R^ is the radius of the cluster. For the moment, we assume R^ to 
be 20/i~^ Mpc, the fall-off point in eq. (0). According to Poisson statistics, the probability 
that this line of sight strikes two or more galaxies in the cluster is 

P(r^) = 1 - [1 + Ar(rx)]e-^('^^). (6) 



Averaging over r±, the probability that a quasar hne of sight strikes two or more galaxies of 
a cluster when it passes through the cluster is J^" P{rj_) ■ 2r±dr±/ R^. Considering that the 
number of clusters the quasar line of sight passes through is 

iVci = / 7TRln,iduil + zf'dz, (7) 

-^min 

where dn = 3000/i~^Mpc is the Hubble radius, and rid is the number density of rich clusters 
today, estimated to be ^ 10^^/i^Mpc"^ (Bahcall 1988). The factor (1 + z)^-^ is the product 
of (1 + zY contributed by the number density of rich clusters at z, (1 + z)~^'^ contributed 
by the differential proper length, and {1 + z)~^ due to the scaling of the proper length at z. 
A universe with the critical density is assumed. The expected number of events for a quasar 
line of sight to strike two or more galaxies in the same cluster, therefore, is 

AfrnuMi = 27r / n,,dn{l + z)''-''dz / P(r^) ■ r^dr^. (8) 

With a typical range of Ziain ~ 1 and z^^^ ~ 3 for CIV absorption lines (SBS), numerical 
integration of eq. yields an event rate of 0.06 per quasar. TVmuiti would remain roughly 
unchanged even if we assume eq. (H) is valid beyond 20 Mpc. It will only increase by 25% 
if we assume that the halo radii of galaxies in the core of clusters are typically 35/i^^kpc 
instead of eq. (H) (while the halo radius of galaxies outside the core remains unchanged to 
maintain the required cross section to intercept quasar lines of sight), resulting a ncore eight 
times higher according to eq. @. Therefore, the expected number of pairs of absorbers with 
200km/sec< AV < lOOOkm/sec due to clusters of galaxies intercepting QSO line of sight is 
< 0.06 per quasar, or 3 for a survey of 55 quasars, which can not explain the ~ 20 pairs 
observed in this velocity range for 55 QSOs (SBS). 
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III. HIGH VELOCITY GASES IN GALAXIES 



It is well known that supernova (SN) shocks or stellar winds can accelerate the interstellar 
medium to velocities as high as ~ lOOOkm/sec. Therefore, the clustered absorption systems 
in the range 200km/sec< AV <600km/sec may originate from the activity of massive stars 
inside the absorber itself, regardless of whether the absorber is a single galaxy halo (Pettini 
et al. 1983), a dwarf galaxy (York et al. 1986) or a mini-absorber inside a galaxy (Sargent 
et al. 1979). A SN typically releases 10^^ ergs in kinetic energy. A 30Mq star also deposits 
~ 10^^ ergs in kinetic energy through mass loss during its main sequence and pulsation life 
time of ~ 5 X 10^ yr (Abbott 1982). Simply on an energetic base, if we assume negligible 
energy loss, the amount of mass a SN shock or stellar wind can accelerate to a velocity of v 
for each SN or massive star is 

lO^^ergs _ , 

In the phase in which SN shocks expand and decelerate to > lO^km/sec, the SN shocks do 
in fact propagate adiabatically with negligible energy loss according to a power law (as given 
by the Sedov solution) (Spitzer 1978) 

V (X (10) 

where t is the time of the propagation. In the case of stellar winds, the gases pushed by 
the strong stellar wind also propagate with negligible energy loss at > 100 km/sec with a 
power law (Weaver et al. 1977) 

V oc r'/^ (11) 

Therefore, for v ~ lOOkm/sec, M„ ~ 10^ Mq. 



We first deal with the case of SN shocks. In our calculations, we assume that gases 
are uniformly distributed in absorbers, and assume that the absorber and the high velocity 
region driven by SN shocks are spherical. The geometry of a high velocity cloud due to a 
SN shock is illustrated in Figure 2. When a quasar line of sight penetrates the cloud, a pair 
of absorption lines are generated with a velocity dispersion 

AV = 2vcose. (12) 

To yield AV > 200km/sec, v has to be larger than lOOkm/sec. The average cross section of 
a high velocity cloud driven by a SN shock to generate a AV > 200 km/sec absorption pair 
before the shock decelerate to lOOkm/sec is 



/■^lookm/scc r /100km/sec\2 

(^)AV>200km/scc = [1 - ( j 



dt 



lOOkm/sec 



vrr 



abs 



tiookm/scc ^ My A^/^Ti /100km/secN^2 



/ x^/'^f-i / iuuKm/sec\ ^-1 dt , , 

yjiZ.) " V ) \ f ' ^^^^ 



^lOOkm/ sec 



where r and v are the radius and the velocity of the high velocity gases respectively. Tabs is 
the radius of the absorber in which the high velocity gases is embedded (which is either a 
single galaxy halo, or a dwarf galaxy, or a mini-absorber inside a galaxy), and Mabs is the 
mass of the absorber, tiookm/sec is the life time of the SN shock when its velocity drops to 
lOOkm/sec. After combining eq. (|]), (|10D , (|12D and eq. (|13D, we get 

I A\ 2 2 /^100km/sec\2/3 

(^)Ay>200km/sec = -vrr^bs( — j . (14) 



and 



, _ 2 Tabs /^100km/sec\l/3 . . 

5100km/secV Mabs ^ ' ^ ' 



where Miookm/sec is aX v = lOOkm/sec. 

Since the survey of SBS showed that ~ 20 pairs of CIV absorption hnes in the range of 
200km/sec< AV <600 km/sec were found in the total of ~ 100 pairs of CIV absorption 
lines (within a sample of 111 CIV absorption line systems), it implies that at any time the 
cross section of these high velocity gases driven by SN shocks to intercept a quasar line of 
light is >10% of that of the absorbers. Therefore, it can be estimated that 

-RsN^ lOOkm/sec (^) Ay>200km/s 



>0.1, (16) 



where i?sN is the SN rate in the absorber. The requirement on i?sN is then 
/lOOkm/secw Mabs \ /lOOkm/secw Mabs \ 



In the case of high velocity regions being driven by strong stellar winds, the only difference 
from the SN case is the power law dependence of v on t, which only results in a modification 
of eq. (|T^) and (|I^) by a factor close to 1. Since massive stars with strong stellar winds 
typically have a mass > lOM©, they inevitably end as SNe. Therefore, our result, eq. (|17|), 
will not be modified significantly. If the cavity created by the stellar wind of a massive star 
is small, e.g., with a radius < 50 pc, the supernova blast of the star later on can catch up 
with the stellar wind and continue to drive the high velocity region according to eq. (|[) until 
radiative losses become important (McCray & Kafatos 1987). This is then an intermediate 
case between a pure stellar wind and a pure SN shock, and the required i?sN will not deviate 
significantly from eq. (|T^ either. 

For single halo absorber models (Pettini et al. 1983), Tabs ~ 70/i^^kpc, Mabs ~ lO^^M©, 
eq. ( [T7| ) gives i?sN > 0.015 per year per lO^^M©, yielding 1.5x10^ SNe in 10^ years (which 
is a reasonable time frame since the detection of these clustered CIV absorption systems 
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ranges from 2; ^ 1.2 to 2; ~ 3 (SBS)). Since each type II SN ejects about O.IMq iron and 
type I SN ejects 6 times as much (Truran & Burkert 1994), this rate translates into at least 
1.5 X 10 ''Mq iron in a 1O"M0 galaxy, or [Fe/H]> -2, in 10^ years. For dwarf galaxy absorber 
models (York et al. 1986), r^bs ~ 8kpc, Mabs ~ 2 x lO'^M©, eq. (0) gives i?sN > 0.0025 per 
year per 2xlO^M0, or 0.12 per year per IO^^Mq, which yields at least 1.2 x IO^Mq iron in 
a 1O^^M0 galaxy, or [Fe/H]> —1, in 10^ years. 

In mini-absorber models (Sargent et al. 1979), the above calculation only applies for 
mini-absorbers with mass Mabs > IO^Mq, since for absorbers comparable to or less than 
10^ Mq, a single SN will disrupt the whole cloud and populations of mini-absorbers will soon 
be destroyed. For mini-absorbers with a larger mass, such as IO^Mq and a radius of 100 
pc, eq. (|T^ yields -Rsn > 10~^ per year per absorber, or i?sN > 5 per year per lO^^M© if we 
consider that 5 x 10^ such mini-absorbers are needed to match the cross section of a single 
70 kpc galactic halo. This SN rate is unacceptably large since the resultant [Fe/H] after 
lO^yr will be larger than 0.6. 

In other words, if the clustered CIV absorption lines in the range of 200km/sec< 
AV <600km/sec are solely due to high velocity gases driven by SN shocks (and/or strong 
stellar winds), the resultant SN rate will be too large for mini-absorber models such as that 
of Sargent et al. (1979), but remain reasonable in terms of the resultant metal yields in the 
single halo absorber model and the dwarf galaxy absorber model. The SN rates calculated 
in the single halo and dwarf galaxy absorber models also lie close to the estimate of the 
current SN rate, 0.005-0.1 per year per galaxy depending on the morphology of the galaxy 
(Cappellaro et al. 1993). (The consideration of different type of supernovae will not change 
our conclusions significantly.) 
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Conversely, with a reasonable SN rate, high velocity gases driven by SN shocks (and/or 
strong stellar winds) play little role in the clustering of absorption lines in the mini-absorber 
models, but may play a significant role in the single halo or dwarf galaxy absorber models. 
Especially when the absorbers had star-bursts at the redshift the absorptions are detected, 
which are expected in the early history of galaxies (Yanny 1990), their SN rates would be 
significantly higher than the SN rate today. The role of SN (and/or stellar wind) driven 
clouds in producing clustered QSO metal absorption lines could then be even more relevant. 

The number of the clustered absorption lines from these SN (and/or stellar wind) driven 
clouds should scale as (i(tiookm/sec(^)Ay>2ookm/sec)/c?'y5 or {AV)~'^, which is much steeper 
than that expected from a galaxy-galaxy correlation function, eq. (|l|). According to this 
dependence law, given the total number of ~ 20 clustered absorption pairs with AV > 200 
km/sec in SBS (1988), the clustered absorption pairs should fall below 1 at AV > 500- 
600 km/sec , which coincides with the observation of SBS (1988) that the clustering of 
absorption lines ends at AV ~ 600 km/sec. Although the low statistics of available data 
and the over-simplification of our present calculations do not allow a realistic comparison 
at the moment, further data with better statistics can potentially test the predicted steep 
dependence of the number of clustered metal absorption lines on the splitting velocity, and 
thereby test the role of SN (and/or stellar wind) driven clouds in producing the clustered 
QSO metal absorption lines. 

Our calculations have been based on a spherical symmetry of the absorber and the 
high velocity region. Although for SNe or stellar winds, a spherical geometry or a nearly 
spherical geometry is most natural, a geometry severely deviates from a spherical symmetry 
(e.g., sheet-like) for these high velocity regions may not be impossible. If so, their cross 
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sections of intercepting a QSO line of sight may increase significantly. However, in such a 
geometry, the velocity is most likely in the direction of stretching of the structures, and the 
velocity projected in the direction of large cross section is small. 

Another concern is the uniform distribution of gases assumed in our calculation. If 
massive stars in absorbers mostly reside in underdensed region, the high velocity regions 
their stellar winds or SN shocks generate will have cross sections larger than that calculated 
from a uniform absorber. Intuitively, however, massive stars should form more likely in 
overdensed regions and gas rich environment. 

IV. SUMMARY 

We showed that while rich clusters of galaxies cannot explain the clustering of CIV 
absorption hues with velocity dispersions in between 200km/sec and ~ 600km/sec, high 
velocity gases driven by SNe and/or strong stellar winds could provide a viable solution to 
these clustered absorptions, provided that the QSO metal-line absorbers are single galactic 
halos or dwarf galaxies. This provides an alternative to the galaxy-galaxy correlation expla- 
nation to these clustered absorbers. Furthermore, the expected number of absorption pairs 
due to high velocity gases driven by SNe and/or strong stellar winds decreases much faster 
with an increasing velocity split AV than that expected from the galaxy-galaxy correlation, 
which provides a testable feature as the statistics of clustered QSO metal absorption lines 
improves. It is also possible that both of these two scenarios, i.e., stellar activities and the 
galaxy-galaxy correlation, work at the same time. Then the significance of the SNe (or 
strong stellar winds) driven gases in providing the clustered absorptions depends on the SN 
rate in absorbers. For single galactic halo models, high velocity gases driven by SN shocks 
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and/or strong stellar winds become important if the SN rate is >0.015 per year per lO^^M©; 
for dwarf galaxy absorber models, these gases become important if the SN rate is ^0.12 
per year per 10^^ Mq. There rates are in or near the range of the estimated SN rate today. 
Therefore, if absorbers had star- bursts at the absorption redshift (between z ~ 1 and ^ ~ 3), 
their SN rates then would be significantly higher than the SN rate today, and the role of high 
velocity gases driven by SN shocks and/or strong stellar winds in providing the clustered 
QSO metal-line absorptions will be important. 
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Figure Captions: 

Figure 1. The geometry of a quasar line of sight passing through a cluster of galaxies. 
Figure 2. The geometry of high velocity gases driven by a supernova shock. 
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